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We investigate the effect of a strong magnetic field on the structure of neutron stars 
in a model with perturbative f{R) gravity. The effect of an interior strong magnetic 
field of about 10^^~^^ G on the equation of state is derived in the context of a 
quantum hadrodynamics (QHD) model. We solve the modified spherically symmetric 
hydrostatic equilibrium equations derived for a gravity model with f{R) = R + aP?. 
Effects of both the finite magnetic field and the modified gravity are detailed for 
various values of the magnetic field and the perturbation parameter a along with 
a discussion of their physical implications. We show that there exists a parameter 



space of the modified gravity and the magnetic field strength, in which even a soft 
equation of state can accommodate a large (> 2 M©) maximum neutron star mass 
through the modified mass-radius relation. 
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I. INTRODUCTION 

Neutron stars, as a remnant of supernova explosions, are an excellent probe of nuclear 
matter in extreme environments. Among the observed neutron-star phenomena, soft 7-ray 
repeaters (SGRs) and anomalous X-ray pulsars (AXPs) are believed to provide evidence 
for magnetars |l|, l2|, ^-c. neutron stars with strong surface magnetic fields of 10^^ ~ 10^^ 
G (see PI for a review). In the interior of these magnetic neutron stars, the magnetic field 
strength could be as high as 10^^ G according to the scalar virial theorem. Such strong 
magnetic fields may affect properties of neutron stars, for example, the relative populations 
of various particles, the equation of state, and the mass-radius relation. Many studies of 



dense nuclear matter in the presence of strong magnetic fields have been reported |4J-|l4|. 
These works have considered the electromagnetic interaction, the Landau quantization 
of charged particles, and the baryon anomalous magnetic moments (AMMs). Indeed the 
detailed analysis of neutron stars with strong magnetic fields is an active area of current 
research. 

Of relevance to the present work are separate studies in which various modifications 
of Einstein's general relativity have been introduced to explain the current accelerating 
cosmic expansion. One can simply introduce a cosmological constant to explain this 
accelerated expansion, but the theoretical prediction for the value of the cosmological 
constant from quantum field theory is many orders of magnitude larger than the value 
inferred from the astronomical data. Equivalently one can postulate negative-pressure 
energy-momentum contribution of unknown origin, i.e. the so called dark energy com- 
ponent with p/p ~ —1, to the matter side of the Einstein equation. An alternative 
to these approaches, however, is to modify the geometry side of the Einstein equation. 
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Such modifications can arise from so called modified gravity theories. A class of modified 
gravity theories is the f{R) gravity (see reviews ISMlSl] and references therein). This 
modified gravity is unique for its simplicity and has been shown to be compatible with 
the constraints from terrestrial laboratory measurements, along with the Solar System 



constraints and neutron star tests 
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II. MODIFIED GRAVITY 



The f{R) gravity theory is simply defined by the following minimal modification to 
the Einstein-Hilbert action, 



S 



167r 



d'^Xy/^f{R) + S'matter, 



(1) 



where g denotes the determinant of the metric g^y and R is the Ricci scalar. Here we set 
G = 1 and c = 1. In the present work we assume that the function f{R) has the following 
linearized perturbative form with respect to a small parameter a without a cosmological 
constant term, 

f{R) = R + ah{R) , (2) 

where h{R) is an arbitrary function of R. The modified TOV equations in this gravity 
are then given by 20l |: 
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where Pa = P + api + •••, and similarly Pa and Mq are expanded in terms of a, i.e. 
Pa = P + aPi + ... and Ma = M + aMi + ... with its zeroth order M = Air J p{r)r'^dr. 
hR{R) is the first derivative of h{R) with respect to i?, whereas h'j^{R) and h'^{R) are the 
first and the second derivatives of hR{R) with respect to r. If we take the a = limit, 
this equation reduces to the standard TOV equation of general relativity. It should be 
noted that we set Ma{M) — )■ 2Ma(2M) compared to the modified TOV equations in [20 1. 



The value of a may be constrained by the observational data. For example, with a 
similar analysis as is done in the present work it has been found in 20| that a < 2x10^ m^. 
For Gravity Probe B it has been deduced that a < 5 x 10^^ m^ 19| , while for the terrestrial 
Eot-Wash experiment the inferred constraint is a < 10^^" m^ 19|]. All of this means that 



allowed values for the parameter a depend heavily on the length scale considered. In 



fact, it can be argued [20| that a length scale given by y/a has nearly the same order of 



magnitude as the typical dimension of the probes used in above tests. Moreover, it is 
also related to the Yukawa correction to the Newtonian potential, ^ exp(— r/A) with a 

Such a violation of the inverse-law of Newtonian 



length scale parameter A = V 6a 
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gravity has been suggested as a possible means to resolve a current dilemma in neutron 



star physics, namely the inconsistency of the super-soft equations of state 23| , obtained 
by interpreting heavy ion collision data FOPI/GSI [24], with the presently observed |25| 
upper mass limit on neutron stars of M > 1.97 ± 0.04 Mq. 

In the present study we consider the combination of strong magnetic fields and modified 
gravity. By way of motivation for this we note that in the five dimensional unification of 
gravity and electromagnetism the Kaluza-Klein action expands into: 

n^ f{R) = i?-a|F|2 , (5) 

where 7?. (R) is the scalar curvature in five (four) dimensions, F is the four dimensional 
electromagnetic field strength and in this case a relates to (the square of) the length 
scale of the extra dimension. In this context, therefore, it is perhaps natural to compare a 
modification of general relativity together with a strong electromagnetic field. Since we do 
not consider charged neutron stars, the Maxwell stress tensor simply reduces to the energy 
density in the magnetic field. Hence, in what follows we consider the effects of magnetic 
field and modified gravity on the neutron star structure combined and contrasted as a 
first step toward a unified picture. As a simplification we consider modified gravity with 
h{R) ~ R^ in Eq. ([2]). Technically, such a term is motivated by Lovelock or Gauss-Bonnet 
gravity in higher dimensions, however, we keep in mind that this could also be loosely 
associated with the Kaluza-Klein electromagnetism. 



III. EQUATIONS OF STATE AND NUMERICAL ANALYSIS 

Many theoretical models for the equation of state (EoS) of the neutron star have been 
developed starting from from an ab initio or effective nucleon-nucleon interaction in order 
to explain the observational mass-radius relation data 126 1 . In this work, we exploit our 

n 

previously developed EoS p^ based upon a QHD Lagrangian and applied to describe the 
interior structure of magnetars. Detailed expressions have been given in that work and 
will not be repeated here. 

Employing this EoS we have numerically integrated Eqs. ([3]) and (jl]). We start from 
the center of the star for a certain value of central pressure, Pc, and then utilize a Runge- 
Kutta scheme with a fixed step size of Ar = 0.001 km. The radius of the star, i?^, 
is identified as the point where pressure drops to a very small value (< 10 dyne/cm^). 
At that point we record the mass of the star, M^. We vary the central density, pc, from 
2 X 10^^ g cm~^ to 1 X 10^® g cm^'^ (to 2 x 10^^ g cm~^ in some cases) in 200 logarithmically 
equal steps to obtain a sequence of equilibrium configurations. We record the mass and 
the radius for each central pressure. This allows us to construct the mass-radius (M-R) 
relation for a given EoS. We have then repeated this procedure for a range of values for 
a to quantify the effect of the perturbative term added to the Lagrangian. 

In the following, we present results for the effect of strong magnetic fields on neutron 
stars with a TOV solution based upon perturbative f{R) gravity. The mass-radius relation 
and the mass vs. central density for each EoS are given for 5 representative values of 
ag = a/lO^cm^ = —2,-1,0,1,2 in the f{R) = R + aR^ gravity. The magnetic field 
strength inside the neutron star is assumed to obey a functional dependence on density 
given by p, yjj 

B{p/po) = B^^^' + Bo[l - exp[(-/3(p/po)1] • (6) 

Here, B^^^^ is the magnetic field at the surface taken as 10^^ G from observations and Bq 
is the saturation value of the interior magnetic field at high densities. In the present work, 
we adopt a somewhat rapidly declining magnetic field strength with density (/3 = 0.02 
and 7 = 3). Since the magnetic field is usually specified in units of the critical field for 
the electron, B^ = ml/e = 4.414 x 10^^ G, the B and the Bq in Eq. (jH]) can be written as 
B* = B/ Bl and Bq = Bq/B^. In this work, we regard the Bq as a free parameter. 

Figure 1 shows the mass-radius relations (left panels) and the mass vs. central density 
(right panels) for the case of a np phase, i.e. no hyperons. The magnetic field strength 



is taken to be B^ ~ 10^~'^ for the 5 representative values of ag. In order to constrain the 
value of a, the recent measurements of the mass and radius of neutron stars, EXO 1745- 



248 m, 4U 1608-52 



and 4U 1820-30 



29| are used. The region bounded by the thin 



)lack line in all M-R plots is the 2a confidence contour based upon these three constraints 
3o|. We also use the mass of PSR J1614-2230 with 1.97 ± 0.04 Mq ^ as a constraint. 



This is shown as the horizontal black line with a grey error-bar. For a viable EoS together 
with an allowed value of the perturbation parameter a, the maximum neutron-star mass 
should lie above this constraining contour. Within the framework of the modified gravity 
considered here these two constraints eliminate many of the possible equations of state. 

The two upper-most figures show results without a magnetic field for the np phase. 
These figures illustrate the effect of the modified gravity alone. Negative (positive) ag 
values give rise to a stiff (soft) EoS. This highlights the interesting feature that modified 
gravity can lead to a mass-radius relation that mimics a stiff or soft EoS for neutrons 
stars. Indeed, if the EoS is ever established to be soft, modified gravity of the sort studied 
here may be required to explain neutron star masses as large as 2 Mq. 

In particular, one may note that negative a values lead to larger masses beyond 2.0 
Mq, while a positive a value tends to diminish the maximum mass. The behavior of the 
EoS in the high density region also has another interesting property in that it shows a very 
strong softness around the region with densities of ~ 2 — 3pc- Magnetic field strengths of 
about Bq ~ 10^~^ have little effect on the maximum mass and mass-radius relations as 
shown in Figure 1. However, with a larger magnetic field of Bq ~ 10^, a stiff EoS and a 
larger maximum mass is obtained as shown in Figure 2. However, if we assume that the 
ctg values are positive, i.e. considering the interpretation of a Yukawa correction to the 
Newtonian gravity, the strong magnetic field effects for Bq ~ 10^ are compensated by the 
modified gravity effects. 

The remaining figure 3 is for the case of a nph phase that includes hyperons. We observe 
exactly the same phenomena in this case, i.e. the effects of modified gravity are seen as 
mimicking a stiff or soft EoS for neutron stars depending upon whether a is negative or 
positive, respectively. In particular we observe that smaller (and more negative) values of 
a allow for a higher maximum neutron star mass. Therefore, the perturbative parameter 
a of modified gravity is a new degree of freedom that can alter the M-R relation and (for 
negative a) causes some equations of state to be viable for neutron star matter which 
might not have otherwise been allowed. 
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IV. CONCLUSIONS 

In summary, in this paper we have considered the effects of strong magnetic fields on 
the neutron star mass-radius relation by using a modified TOV equation derived from the 
f{R) = R + OiR^ modification of Einstein's general relativity. We fixed the upper bound 
of the expansion parameter ag = a/lO^ cm^ of the gravity model using several constraints 
based upon neutron star observations. It turns out that in the case of ag > 0, effects by 
the modified gravity can be compensated by those of a strong magnetic field. However, 
the case of ag < is completely different: in this regime some equations of state which 



were not viable for neutron stars in the case of genera 



ativity become viable again. 
This is in accord with the previous results given in |20l . l3l| on the consequences of the 
perturbative modification of general relativity. 

It would be instructive to repeat the present analysis for the case of a non-spherically 
symmetric, but axially symmetric space-time, because a high magnetic field could alter 
the almost spherically symmetric shape of a neutron star to an axially symmetric one. 
Then one would better understand the interplay between the high magnetic fields and 
the effects of gravity on the physics of neutron stars. In the case of white dwarfs with 
high magnetic fields, the magnetic field's possible effect on the star's geometry and the 



consequence for its mass and radius has been commented upon in 32]. The present work 
should be considered as a first step toward that line of research. 
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FIG. 1: (Color online) Mass-Radius relation (left panels) and the mass vs. central density 
relations (right panels) corresponding to an EoS with a np phase and magnetic fields with 
Bq = 10^~'^ in the modified gravity model characterized by values for the parameter a. Here 
"EoS = nbm b-np" stands for the EoS with Bq = n x 10*" in a np phase. The observational 
constraint [30] on the mass-radius relation is shown by the thin dashed black contours. The 



mass of M = 1.97 ± O.O4M0 for PSR J1614-2230 
grey error bars. 
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FIG. 2: (Color online) Same as Fig.l but with magnetic fields Bq = {2 ^ 5) x 10^ 
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FIG. 3: (Color online) Same as Fig. 2 but corresponding to an EoS with a n-ph phase, i.e. 
including hyperons. 
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